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We elucidate the topological features of the entanglement entropy of a region in two dimensional 
quantum systems in a topological phase with a finite correlation length ^. Firstly, we suggest that 
simpler reduced quantities, related to the von Neumann entropy, could be defined to compute the 
topological entropy. We use our methods to compute the entanglement entropy for the ground 
state wave function of a quantum eight-vertex model in its topological phase, and show that a finite 
correlation length adds corrections of the same order as the topological entropy which come from 
sharp features of the boundary of the region under study. We also calculate the topological entropy 
for the ground state of the quantum dimer model on a triangular lattice by using a mapping to 
a loop model. The topological entropy of the state is determined by loop configurations with a 
non-trivial winding number around the region under study. Finally, we consider extensions of the 
Kitaev wave function, which incorporate the effects of electric and magnetic charge fluctuations, 
and use it to investigate the stability of the topological phase by calculating the topological entropy. 



I. INTRODUCTION 

Different phases of condensed matter are usually dis- 
tinguished by broken symmetries of the physical system 
measured by local order parameters, an approach that 
is ultimately justified by the renormalization group. In 
quantum systems at zero temperature, these local or- 
der parameters correspond to local observables having 
nonzero expectation value in the ground state wave func- 
tion in the broken symmetry phase. However, quantum 
systems can also display phases where the character of 
the state is not measurable by any local order parame- 
ter. An example is topological order— i^, a form of quan- 
tum order marked by a sensitivity of the ground state 
degeneracy to the topology of the system in the absence 
of any broken symmetry. An experimental system with 
topological order is the two-dimensional electron gas in 
the regime of the fractional quantum Hall effect li A num- 
ber of simple and very idealized theoretical models^ii^i^ 
(not yet realized in experiment) with topological phases 
have now been constructed. 

Recently, Kitaev and Preskilfi and Levin and Wen^ 
proposed that, in two-dimensional systems, whether a 
ground state wave function j^I') has topological order 
can be determined by computing a nonlocal quantity 
called the topological entanglement entropy —7. The 
computation involves dividing the system into two sub- 
regions, A and B; obtaining the reduced density matrix 
Pa = T7-b|5')(\1/| by tracing over the degrees of freedom 
in B; and then calculating the von Neumann entropy 
Sa = —TtaPa In PA by tracing over region A. The ob- 
servation was that if the boundary between regions A 
and B has perimeter L and is sufficiently smooth, then 
for a system with short range correlations, the entropy 
should have the scaling form: 

SA=-aL--i + ... (1.1) 

where a is a non-universal finite coefficient. Here ". . . " 



represent terms that vanish in the limit of a large re- 
gion, L — > (X). The subleading constant term, the topo- 
logical entanglement entropy^ (or just the "topological 
entropy"), is universal: it depends only on the type of 
topological order. More specifically, they find^i^ 

7 = In 2? (1.2) 

where 2? is the so-called total quantum dimension. For 
a normal state, 2? = 1 while for topological orders de- 
scribed by discrete gauge theories, V is the number of 
elements in the gauge group. In the general case, V need 
not be an integer but in this paper we will concentrate on 
the simplest type of topological order, namely Z2 gauge 
theory where V = 2. Examples of models with Z2 topo- 
logical phases include the Z2 gauge theory and its ex- 
treme deconfined limit, Kitaev's toric codo^, the quan- 
tum eight vertex model^, and the quantum dimer model 
on the triangular— and other latticesjiS— 

Eqs. (|l.ip and (|1.2p were derived in Rcf. 7i using topo- 
logical quantum field theory, believed to capture the uni- 
versal properties of topological phases, while in Ref. H, 
these expressions were obtained for the collection of 
"string- net condensate" wave functions, scale-invariant 
wave functions, believed to be "fixed point" models of 
topological order. A common feature of both approaches 
is that the correlation length is (implicitly in Ref. 0) ex- 
actly zero. 

In the present work, we study systems with finite cor- 
relation length ^ to understand the effect of local correla- 
tions on the subleading part of Eq. (jl.ip and the issues in- 
volved with extracting the topological entropy from non- 
universal corrections. Our motivation is twofold. The 
first is to understand the sense in which 7 is a univer- 
sal property of a topological phase. The calculations in 
Refs. andU are based on the assumption that for sys- 
tems with a finite correlation length ^, the universal be- 
havior should apply on length scales much longer than ^. 
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In this paper we reexamine the effects of the existence of 
a finite correlation length in the context of specific mod- 
els whose wave functions are known for all ^, and verify 
that, indeed, this behavior is a property of the topological 
phase and not just of the ^ ^ limit. 

Having a finite correlation length should change the 
value of the non- universal coefficient a in Eq. p.ip as well 
as give non-universal corrections to the constant term. 
However, Refs. and H offered strong additional argu- 
ments for why it should still be possible to define a uni- 
versal topological entropy. In Ref. 0, it was noted that 
correlation length related corrections to the subleading 
term would come from parts of the boundary that were 
not smooth (on the scale of ^) such as corners. This led 
them to calculate the topological entropy by partition- 
ing the system into subregions and considering a special 
sum of von Neumann entropies designed to exactly cancel 
contributions related to the boundary and corners leav- 
ing only the universal part 7. In Ref. it was noted 
that a defining feature of topological order is nonlocal 
correlation revealed by the existence of closed loop op- 
erators having nonzero expectation value. They noted 
that having a finite correlation length would amount to 
a redefinition of the relevant loop operators, including 
perhaps an effective "fattening" of the loop. Either way, 
the von Neumann entropy of an annular region would be 
lower than the value expected for purely local correla- 
tions and the difference approaches a universal value (27 
for their geometry) in the large L limit. In this paper, we 
will show how these ideas work by explicit calculations 
on models with finite correlation length. 

Our second motivation is that the practical importance 
of the topological entropy as a diagnostic for topological 
order is limited by the ease with which the universal part 
can be separated from finite size effects and other non- 
universal corrections due to finite correlation length, and 
(on a lattice) ambiguities in defining the length of the 
boundary. These latter two issues make the process more 
subtle than just calculating Sa versus L and then find- 
ing the intercept. While such effects can be accounted 
for by the constructions mentioned abovejii^ there is still 
the problem of attaining a large enough system size that 
the scaling form Eq. (|l.ip becomes accurate. In addition, 
if the calculation involves approximate methods, then a 
natural issue with these constructions is that 7 is a sub- 
leading term obtained by adding and subtracting large 
numbers. This can make approximate calculations diffi- 
cult as the errors can be comparable to the quantity of 
interest. 

These issues were confronted in recent numerical works 
aimed at verifying Eq. (jl.2p for the simplest examples 
of topological phases with finite correlation lengths. In 
Refs. , 12il3i . 7 = In 2 was verified for the toric code^ and 
triangular lattice quantum dimer model at the Rokhsar- 
Kivclson (RK) point^ii^ respectively, by directly com- 
puting von Neumann entropies from the exact reduced 
density matrix. Similar calculations, performed for frac- 
tional quantum hall wave function a^^i^^ , provided values 



of 7 consistent with Eq. (|1.2p . but with errors due to 
finite size effects. In this paper, through model calcu- 
lations, we aim to clarify how topological information is 
"stored" in the ground state wave function and how a 
finite correlation length affects the entanglement. This 
understanding will lead us to propose reduced quantities 
which should contain the same topological information as 
the von Neumann entropy but may be simpler to work 
with numerically. While our calculations will apply for a 
particular class of wave functions, our conclusions should 
hold more generally. 

In this paper, we mainly concentrate on wave functions 
whose amplitudes have the same form as local Boltz- 
mann weights of related models in 2D classical statis- 
tical mechanics.— i^i^iii^'^ In Section |llj we will show that 
for such states, calculating the von Neumann entropy of 
region A is equivalent to determining probabilities for 
configurations of the boundary of A. (A closely related 
problem, the entanglement entropy at z = 2 conformal 
quantum critical points in 2D, was considered in Ref. [isl . 
which the present work generalizes.) From this perspec- 
tive, in Section IIIII we will consider the von Neumann 
entropy of the ground state of the quantum eight-vertex 
model near the Kitaev point^, where the ground state 
wave function is equivalent to that of the toric code^. 
The Kitaev point is located inside a phase with Z2 topo- 
logical order and at this special point, the correlation 
length is zero. Away from this special point, the corre- 
lation length increases in a known manner. In perturba- 
tion theory, we verify the assertion of Ref. that finite 
correlation length corrections to Eq. (jl.ip come from ap- 
propriately defined corners of the boundary of the region 
under study. In Section ITvl we consider the ground state 
of the triangular lattice quantum dimer model at the RK 
point^iii. We show that the topological entropy can be 
related to the emergence of topological winding sectors 
in the calculation of the entanglement entropy. Finally, 
in Section |V] we consider the effects of topological defects 
on the entanglement entropy of the quantum eight ver- 
tex model. We first show that mobile defects in the wave 
function, which drive the system into a confining phase, 
make the topological entropy vanish. However, in the 
perturbative regime where defects are slightly favored, 
topological entropy is shown to remain robust in low or- 
ders in the perturbation expansion series, and only cor- 
rections due to the boundary shape are allowed. Along 
these lines, we identify the terms in the wave function 
which add topological corrections to the topological en- 
tropy, and when proliferating, should lead to a topolog- 
ical phase transition. Our results clearly suggest that, 
within the radius of convergence of perturbation theory 
(which is indeed finite) , and for regions much larger than 
the correlation length ^, the topological entropy 7 is a 
universal (and hence constant) property of the topologi- 
cal phase. In Section IVll we conclude by discussing open 
problems and implications of these ideas for numerical 
calculations of the topological entropy. Finally, in Ap- 
pendix |X] we present a summary of known results of the 
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"Li gauge theory relevant to this discussion. 

II. ENTANGLEMENT ENTROPY AND 
MEASURES FOR TOPOLOGICAL ORDER 

In this Section, we discuss the problem of calculating 
the von Neumann (entanglement) entropy for a class of 
wave functions that includes the ground states of various 
two-dimensional models of the type considered in this pa- 
per. The key property of these wave functions is that the 
amplitudes of the various states are local functions of the 
degrees of freedom, i.e. the amplitudes have the form 
of a Gibbs weight in a suitably defined classical statisti- 
cal mechanical system. Consequently, the norm of such 
wave functions has the form of a partition function of 
the related classical statistical mechanical system. The 
main aim of this Section is to show that the calculation 
of the entanglement entropy in such wave functions has 
an analog in the equivalent classical statistical mechani- 
cal problem. This strategy, applied to generic topological 
phases with finite correlation length, leads to natural gen- 
eralizations of the relations recently derived in Refs. [isl 
and[li. 

We begin by dividing the system into simply connected 
subrcgions A and B. Each degree of freedom in the 
system belongs to either region A ox B and the closed 
boundary curve separating these regions is denoted by 
r. A convenient basis to work with is the set of product 
states of the form |a, 6) = |a) ® where |a) and |6) de- 
note particular configurations of the degrees of freedom 
in regions A and B respectively. This basis is complete 
and orthonormal. 

In this basis, an arbitrary wave function will have the 
form |\1>) = ^ T,a,b '^a,b \a, b), where the Ca,bS are com- 
plex amplitudes, and Z — J2a t ka^fcP ^^e normaliza- 
tion of the state l^"). 

The wave functions we are interested in have three ad- 
ditional properties: 

1. The sum is over a restricted set of states that may 
be written in terms of the action of a set of transfor- 
mations G on some reference state |a, b) = g |ao, 60)1 
where g G G- 

2. For the states allowed in the sum, we assume the 
probability amplitudes are local, that is: 

Cab = Cah (2.1) 

i.e., a change of the state in region A is directly 
connected to the corresponding change of the am- 
plitude, irrespective of the state in region B. 

3. We assume the set of transformations G form a 
group, i.e., The elements of G satisfy the four group 
axioms: (a) closure. If (71,52 S G, then the com- 
posite transformations g^ = gig2 and (73 = 5251 arc 
also in G (though 33 need not equal 53). (b) as- 
sociativity. For all gi, 52, and in G, {9192)93 = 



91(9293) ■ (c) the identity transformation is in G and 
(d) If 5 e CJ, then so is the inverse transformation 

These properties imply that our wave functions of inter- 
est have the form: jvl/) = SgGg ^a'^b9 I'^o, ^0), where 
the reference state |ao, &o) and the group G encode what- 
ever constraints the wave function respects. The nature 
of the group of transformations G depends not only on 
the system but also on the representation chosen to de- 
scribe the states. In the case of the systems that we 
discuss in this paper, the quantum dimer model, the Z2 
gauge theory and the quantum eight vertex model, there 
is a local Z2 gauge symmetry. In all these systems there 
are two mutually dual representations, both of which as- 
sociated with the group Z2, although the physical inter- 
pretation will be different. As shown in Appendix lAl in 
the electric representation, the Hilbert space is the space 
of a set of loop( "strands") configurations. In this case, 
the Z2 group of transformations {g} is the Braid group 
with the associated Temperley-Lieb algebra, represented 
in the Hamiltonian by the flip (or resonance) term. In 
the language of the Z2 gauge theory, there is a one-to-one 
correspondence between these transformations and the 
set of Wilson loops of the gauge theory, Wj = Yieer '^i' 
where Tj arc closed loops and {£} arc the links on that 
loop. Conversely, in the magnetic representation, the Z2 
group is just the group of local gauge transformations. 
In this picture the gauge transformations are in one-to- 
one correspondence with the set of dual Wilson loops 
Wk = Yli^f^ ' "^here F^ are closed loops on the dual 
lattice, and {£} are links of the direct lattice intersected 
by Ffe . It is a peculiarity of systems with a Z2 symmetry 
that these dual representations are isomorphic to each 
other. The electric and magnetic Wilson loops obey a 
non-trivial algebrai^ 

For example, in the quantum dimcr model, G is the set 
of operators whose action is to move dimers from occu- 
pied to empty links along flippable loops. In other words, 
G is the set of operations which allow one to transform a 
given dimer covering of the lattice into another without 
ever violating the constraint that a site may have one 
and only one dimer. Similarly, in the quantum eight ver- 
tex model, G is the set of operators whose action is to 
reverse the directions of arrows along flippable loops. In 
other words, G is the set of operations which allow one 
to transform a given eight-vertex covering of the lattice 
into another without violating the constraint that a site 
can only have 0, 2, or 4 arrows flowing into it. 

We define the reduced density matrix of region A to 

be 

PA = TrB\^){^\ (2.2) 
and its von Neumann entropy is defined as: 

Sa = -Tr^ [pA log pa) (2.3) 
A useful way to calculate Sa is through the so-called 
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"replica trick"^i2a: 



Sa^- lim — Ttap'X = lim (2.4) 

n~*l dn \ J n^l 

where p\ is the n-th power of the reduced density matrix 
PA- 

It is an implicit assumption of this "replica" approach 
that there is a well defined and unique analytic contin- 
uation of the generalized "entropies" SA,n defined for n 
in the set of positive integers to a region of the complex 
plane which includes n ~ 1. As wc shall see, this assump- 
tion holds for the cases we are interested in. (We note, 
however, that there are systems, studied in the recent 
literature fS^i^ for which the existence of the n ^ 1 limit 
is quite non trivial.) 

For a general wave function, TtaPa is given by: 

Z^'Ti-APl = 

^ai,biCa2,bi^a2M^a3,b2 ' ' ' ^an,bn^ai,bn 

ai,bi ....,a-n, ,bn 

(2.5) 

For a constrained wave function, c* i will be zero un- 

' (12,01 

less |a2,6i) is one of the allowed configurations in the 
ground state sum and so on. If in addition, we assume 
the condition (|2.ip . then the expression may be written: 



T^^APl 



- E 

a\ ,bi---arL -b-n 



\CaArCtA'\Ca2\' ■ 



\Ca,f\hS 



(2.6) 



where the prime denotes that the sum over bi is restricted 
to states where both |ai, bi) and |a2, bi) are allowed con- 
figurations and likewise for 62 and so on. If we view the 
labels l,2...,n as denoting n different copies of our two- 
dimensional system, then the sum in Eq. (|2.5p can be 
interpreted as a partition function for n systems where 
region A of copy -I- 1 is constrained to be consistent 
with region B of copy k and so on. 

An additional simplification arises if the allowed states 
in the wave function are related by the group property. In 
this case, we have "transitivity", i.e. if |ai,6i), |a2,6i), 
and |a2, ^12) are allowed configurations, then so is |ai, 62)- 
The implication is that the sum in Eq. (j2.6p is the parti- 
tion function for n copies that are constrained such that 
region A of any copy is consistent with region B of any 
other copy. We may write this as: 



T^APl 



Z[n] 

z[iY' 



(2.7) 



where Z[n\ is the partition function of n copies con- 
strained to agree with each other on the boundary F (in 
the sense discussed above) but are otherwise independent 
and Z[l] = Z is the partition function of a single copy. 
This result was derived in Ref. [l^ by a different line of 
arguments. 



The most trivial way that the n copies agree on the 
boundary F is if the state of the boundary is the same in 
all of them. However, this is not the most general way 
of agreement. In general, as we noted above, the opera- 
tors g can be regarded as either products of Wilson loop 
operators Wj, or of the dual Wilson loops, Wk, of the 
corresponding gauge theory. Such a Wilson loop opera- 
tor might involve degrees of freedom in region A (PF^"^'), 
m region B (VF(-^)), or both (IF(^^)). What is clear now 
is that the agreement of all n copies persists if no VF^"*^) 
operators transform any of the copies. This means that 
the degrees of freedom on the boundary F might change 
through application of VF^"*) or IF^^) and the agreement 
of the n copies is not affected. The loop operators VF^"^^) 
naturally form equivalence classes, each defined by the 
set of loops W^^^^ modulo contractible loops defined in 
either region A or B but not in both. Each class is la- 
belled by an irreducible loop W^^^\ However, if any 
of the possible irreducible operators VF^^^) is applied on 
any of the copies, then the agreement is violated and it is 
restored only if the same IF'^^'^^ or any of its reducibles 
is applied on all the other copies. 

As a concrete example, consider the quantum dimer 
model where the boundary loop F is drawn on the direct 
lattice (see Fig. |3(b)] ). We consider the lattice links inside 
and along the boundary as belonging to region A and 
the links outside as part of region B. As will be shown 
in Section IIV( we can classify dimer configurations by 
looking at the lattice sites which lie along F and for each 
of these sites, specifying whether its dimer belongs to 
region A or B. Two dimer coverings that agree along F 
in this mannet^ may be said to "agree on the boundary" 
in the sense discussed above: if attempt to make a new 
dimer covering by combining the dimer pattern of region 
A of the first copy with that of region B of the second 
copy, the new configuration will automatically satisfy the 
hardcore dimer constraint on the boundary. It is easy to 
see that flipping the dimers from occupied to empty links 
along loops which live solely on links of region A ot: B^ 
i.e. acting on the configuration with operators IF^"^) or 
W^^'^ in the above paragraph, will not change the state of 
the boundary as just defined. However, performing this 
operation along a loop which intersects the boundary, i.e. 
acting on the configuration with operators IF^t-^^' in the 
above paragraph, will change the state of the boundary. 

Returning to the general case, it is useful to classify 
configurations based on the states of their boundaries. 
This procedure leads to the following identification: 



E 

se{G(-4B)} 



mi 



(2.8) 



To explain the notation, we choose a reference config- 
uration which also defines a reference configuration of 
the boundary. As mentioned above, we can then group 
the various transformations G*-^^^ which will change the 
boundary configuration into equivalence classes based on 
the final state of the boundary. We denote the collection 
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of equivalence classes by {G^"^^^} and within each class, 
we may select some 5 as a representative element. The 
sum in Eq. 12.81 is over the set of these representative g's, 
which is an irreducible set. The partition function ^[l]p 
is defined as 

(2.9) 

where a, b run along all possible configurations of any of 
the copies which are produced by applying products of 
operators, respectively, on a reference con- 
figuration. In other words, ^[l]p is the partition function 
of a single copy of the system (note Z[l] as opposed to 
Z[n\) where the configuration of the boundary T is now 
fixed to the configuration obtained when transformation 
g acts on the reference state. 

In this way, if wc use (|2^ . ([ZT]) and (fZS]) . we have: 

Sa^- E PflogW) (2.10) 

where = Z[l]f/Z[l]. 

From another point of view, given that the sum in Eq. 
(|2.10[) runs over all possible irreducible (either gauge or 
Temper ley-Lieb, depending on the case) transformations 
g & {G'"^^-'} of the boundary configuration, involving 
degrees of freedom in both regions A and B, we have: 

^--('-(ff)). 

The form of Eq. [2TT] is quite instructive, showing that the 
entanglement entropy is an average quantity over gauge 
transformations, as one should expect, since the bound- 
ary has no special properties. 

In a topologically ordered state (non-critical) which is 
defined by a ground state wave function which satisfies 
the property in Eq. (|2.ip . the calculation of the entangle- 
ment entropy leads to the general resultii^,: 

SA = a{0L + b{0~^ + O{l/L) (2.12) 

where a(^),6(^) are correlation length dependent con- 
stants. By the argument of Kitaev and Preskilli and 
Levin and Wen^, the constant term 7 is expected to 
be a topological invariant and as such it should remain 
constant throughout the topological phase. Its value 
is known in the topological limit , where the correlation 
length vanishes, to ho^^ 7 = In where di is the 

quantum dimension of particles of type i and the sum is 
over all non-trivial supcrselection sectors of the medium. 

On the other hand, the average over all possible trans- 
formations in Eq. (|2.1ip shows that "measuring" the en- 
tanglement entropy of a region of the medium does not 
endow its boundary with any special property that would 
make it different from any other place in the system. Mo- 
tivated by these observations, we note that since 7 is a 



topological invariant, its value should in fact be the same 
for all members of the ensemble of transformations im- 
plied in Eq. f2.1ip . Thus, it should be possible to obtain 
this quantity by "fixing completely the gauge" on the 
boundary between regions A and B which should not af- 
fect the topological properties of the state. In this way, 
we define the following reduced quantity: 

where go is a single choice of all possible G'-^^-'s. We 

(s) 

conjecture that S\ contains the same topological infor- 
mation as 5*^. In other words: 

S^;^ = ~a{OL + biO-l + d{l/L) (2.14) 

The calculation of this quantity involves an effective split- 
ting of the system, similar to the procedure that was used 
in Ref. H^. On the other hand, the fact that S'^^'' has 
the form of Eq. (|2.14p can be trivially proven for topo- 
logical systems with zero correlation length, such as the 
ones m Refs. Moreover, we verify in the follow- 

ing sections that such a form still holds for systems with 
a finite-correlation length, such as the quantum eight- 
vertex model of Ref. in its topological phase, in the 
deconfined phase of the Z2 gauge theory, and in topolog- 
ically ordered ground states of quantum dimcr models at 
the RK pointsiiii. 

From a numerical point of view, the advantage of this 
reduced quantity is that the whole set of matrix elements 
of the density matrix of the full system is not necessary. 
Instead, one needs the density matrix of a specially local 
wave function, making the task of extracting 7 from a 
wave function much simpler. Moreover, the constructions 
by Levin- Wen and Kitaev-Preskill can be used for S)^ 
as well, but now the number of additions-subtractions of 
large numbers that are needed is minimized. 

The use of the representation of entanglement entropy 
of Eq. (|2.11|) and the consequent reduction of Eq. ()2.13p 
is formally possible only when Eq. (|2.ip holds. Eq. (j2.1|) 
expresses the requirement that regions A and B arc inter- 
acting ultra-locally, only through the degrees of freedom 
of the boundary F. However, more generally, one should 
consider ground state wave functions where the ampli- 
tudes Co, 6 are defined as Ca,b = CaCf,cJj"*(d), where c"J,'(d) 
depends on the relative configurations of regions A, B 
at a distance d from the boundary F. The distance d 
characterizes microscopic properties of the ground state 
wave function and typically d > a where a is the lattice 
spacing. The distance d should not be confused with the 
correlation length of the system ^ which could become in- 
finite. Instead, it should be viewed as a consequence of lo- 
cal short-range (of the order of d) interacting terms of the 
Hamiltonian, whose ground state we consider. To make 
Eqns. (|2.11|) and (|2.13p applicable in these cases, the def- 
inition of the reduced density matrix should be modified. 
If |F*) represents a generic configurational state of the 
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region T near the boundary (at distance d) in region A, 
then we define: 



(r,|TrBp|r, 



(2.15) 



In this way, all degrees of freedom in region B are in- 
tegrated out and the degrees of freedom in region V are 
fixed, p"^ is a reduced density matrix of region A. Now, 
the generalized entanglement entropy of region A using 
as the reduced density matrix is : 



Pa ^"^Pa 



(2.16) 



where the sum runs along all possible configurations of 
region F. Now, Eqns. (|2.7p and ()2.8p can be applied with- 
out any change, because the property in Eq. (|2.ip holds 
for regions A, B and thus, a partition function represen- 
tation is possible. Following the same type of reasoning, 
we end up to the following relation: 



^A 



In 



^i(r) ^i(r) ' 

A B 



m 



(2.17) 



where now i(F) labels configurational states of region F. 
Following the same type of reasoning as above, we argue 
that, if io(F) labels a specific configuration of region F, 
then the reduced quantity: 



^A 



In 



A B 



Z[l] 



(2.18) 



contains the same topological information as Sa- 

Even though such a generalized quantity, where the 
boundary is effectively "fattened" , should be used for the 
case of general wave functions, in the case of the quantum 
eight-vertex model or quantum dimer models at the RK 
point this is not necessary. 



III. QUANTUM EIGHT- VERTEX MODEL 
NEAR THE KITAEV POINT 

In this Section we compute the entanglement entropy 
of the ground state wave function of the quantum eight- 
vertex model of Ref. in its quantum disordered, topo- 
logical, phase. We will follow the methods discussed in 
the previous Section. 

The quantum eight- vertex model of Ref. is a a two- 
dimensional quantum generalization of the classical two- 
dimensional eight vertex model^^, which is integrable. In 
the classical model, the degrees of freedom are arrows 
that live on the links of a square lattice; arrows on hori- 
zontal (vertical) links point left or right (up or down) . In 
addition, there is a local constraint that each vertex may 
have only an even number of ingoing (or outgoing) ar- 
rows. This results in eight types of allowed vertices, (cf. 
Fig 1(a) I The partition function of the classical model is: 




FIG. 1: (a) Allowed vertices in the classical eight- vertex 
model, (b) The phase diagram of the classical eight-vertex 
model. The phases I and II represent "antiferroelectric " 
phases dominated by d and c vertices respectively. The mid- 
dle D phase is disordered. The point K is where c = d = 1. 
In the quantum version of the model, the ground state wave 
function at this "Kitaev point" is the same as the ground 
state of Kitaev's toric code^. 



(3.1) 



Due to global arrow neutrality (the ingoing equal the 
outgoing arrows), the weights a, b are redundant, and can 
be set to I. The phase diagram of this model in terms of 



the weights c and d is shown in Fig 1(b) There are two 
ordered phases for > cP + 2 and (P~>c^ + 2 separated 
from a disordered phase by lines of continuously varying 
exponents (c^ ~ (P + 2 and (P = + 2). Moreover, for 
c = 0, < 2 and d = 0, < 2, there are additional lines 
of fixed points. The region of interest for our calculation 
is the vicinity of the point c = d = 1. 

The quantum generalization of this model begins by 
placing Pauli matrices cr^ on the links of the square lat- 
tice where k labels the link and a = 1,2,3 labels the 
Pauli matrix. We work in the representation where cr^ 
is diagonal and up spins correspond to arrows pointing 
up or right. In this language, the eight-vertex constraint 
corresponds to: 



-x/2'^l+x/2'^l+y/2'^l-y/2 1^) ~ 1^) 



(3.2) 



where x and y are lattice unit vectors. This expression 
holds for any vertex r and |G) is a state within the mani- 
fold spanned by the eight vertex configurations. We then 
consider the following Hamiltonian: 



^ 5^ '^r+x/2'^r+x+y/2'^r+x/2+y'^r+y/2 (3-3) 



where the sum is over plaquettes. The individual terms 
are operators that fiip the spins (or equivalently the ar- 
rows) around a plaquette. These flip operators commute 
with the constraint so it is possible to diagonalize both 
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simultaneously. As shown in Refs. and[^, for c = d = 
this results in a ground-state wave function given by: 



\G 



K 



(3.4) 



{C} 



where C corresponds to an eight-vertex configuration of 
the square lattice. This wave function is exactly the 
ground-state wave function of Kitaev's toric code. The 
norm of this wave function corresponds to the partition 
function of the eight- vertex model at c = d = 1, which 
we label the Kitaev point (K in Fig. 1(b) ). Moreover, in 
Refi, it was shown that the Hamiltonian Hk can be lo- 
cally deformed in such a way that the ground state wave 
function now becomes: 



1 



(3.5) 



{C} 



where now c, d are free parameters of the Hamiltonian 
and: 



^r~x/2 



(3.6) 



above, and explained in Appendix [XI the flip operator is 
a representation of the Temperley-Lieb generator acting 
on the loops (or strands) of the eight-vertex model. 

Now, we can label each plaquette with an Ising vari- 
able Tj , living at the center of each elementary plaque- 
tte, which will indicate in an eight-vertex configuration 
whether the plaquette- flip operator is applied or not. In 
the vacuum state |0), Tj = — 1 for any plaquette. In this 
representation, we have: 



Iff 

Iff 
'^T+x/2 



Q Q 
T- ~T~ 
r—x r 



'r,3 'r,4 



(3.10) 
(3.11) 



where the indices 1, 2, 3, 4 are defined in Fig 2(b) In this 
representation, the density variables nc,nd p.6p . (|3.7ll . us- 
ing the definitions (jS.lOp . (|3.1ip can be expressed as fol- 
lows: 



n-dig) 



1 

4^ 



l-l I T-ff _ _ff _ _s _ff _ff _ff 



'''r,2''r',3 + '''r, 1 '''r,4 '''r ,l'''r ,2'''r ,3'''r ,4) 



(3.12) 



fidiC) 



16 ^ 



x/2 



- ^r+ji/2 - ^r- 



y/2 



(3.7) 



The norm of this state, (Gsy IGgy), is just the parti- 
tion function of the classical two-dimensional eight-vertex 
model, Eq. p.ip and thus, the ground-state phase dia- 
gram is identical to the classical one. The main differ- 
ence is that what in the classical system is the disordered 
phase, in the quantum system is the topological phase. 

We also may rewrite the ground-state wave function 
p.5p in a more helpful and interesting way: 



1 



8V) 



and its norm: 



:^c"'=(s)d"''(9)5|0) 



{ff} 



{g} 



(3.8) 



(3.9) 



where |0) represents the "vacuum" state, where |0) = 
|0) for any r and \g) = g\0) corresponds to the repeated 
actions of the flip operator on the vacuum state. More 
specifically, we may think of 5 as a product of single 
plaquette-flip operators Y[a "'r ^^'^ then, \g) differs from 
|0) only in that some plaquettes have been flipped. It 
is important to notice that the single plaquette-flip op- 
erator can be applied to all but one plaquette, because 
the product of all plaquette-flip operators on a closed 
surface is equal to the identity operator. As we noted 




(a) 



(b) 



FIG. 2: (a) A pictorial representation of the boundary line 
r and the dual lattice variables which have to remain frozen 
in the calculation of the entanglement entropy, (b) The dual 
lattice variables with respect to a lattice site of the direct 
lattice. 

These expressions will simplify the calculation of the 
entanglement entropy. 



A. Entanglement entropy near the Kitaev point 

To compute the entanglement entropy in the disor- 
dered phase of the eight-vertex wave function we begin by 
dividing the eight-vertex lattice into two regions A and 
B, separated by a single closed boundary line F on the 
dual lattice (Fig. 



2(a)). We will follow the approach of 



Sectionini The degrees of freedom are arrows that live on 
the links in region A or B. As shown in Fig. 



2(a) 



region 

A contains all the degrees of freedom (links) that define 
the boundary, but the r variables living at the dual sites 
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around the boundary, pointed out in the figure, have to 
be fixed in the calculation of pp, according to the defini- 
tion we gave for pp. An important feature of this choice 
of the boundary is that the species {i.e. whether or not it 
is a c vertex) of a vertex in region A is determined solely 
by the links in region A and similarly for the vertices of 
region B (the eight vertex constraint means that knowing 
three links automatically determines the fourth). There- 
fore, the probability weights of the wave function have 
the property (j2.ip which was needed in order to think of 
the calculation in terms of boundary probabilities. 

At the K point (c = c? = 1), for any choice of the r 
variables on the boundary, the probability pp is: 



q 

Pr 



Nr-1 



(3.13) 



where Ns is the total number of sites in the (direct or 
dual) lattice; A'p is the number of direct lattice links or 
dual lattice sites on the boundary, and the —1 originates 
in the global constraint: 



1 



(3.14) 



This result leads to the well known result for the en- 
tanglement entropy of the Kitaev's wave function: 

5*^ = -^Ppnpj^ = (In2)i-ln2 (3.15) 



As soon as c 7^ 1 and d 1, the correlation length 
^ of the ground-state wave function becomes non-zero. 
In particular, as the critical line is approached on the 
left or right, the correlation length is known^l to diverge 

as e ||c - d| - 2r^/|''^*''""'(^)|. We are interested 
in the effects that the correlation length has on the en- 
tanglement entropy. Now, each transformation g has a 
"Gibbs weight" which makes it more or less favorable 
than others. We argue that any other way of introduc- 
ing a correlation length is going to cause qualitatively 
similar effects to the ones that we describe here. The 
reason is that a correlation length is expressed gener- 
ally by the fact that some eigenstates in some basis are 
more favored than others. We can study these effects in 
what in terms of the related classical eight-vertex model 
amounts to the (rapidly convergent) high-temperature- 
type expansion (we will consider the case c l,d = 1). 
We have: 



cosh(i Inc^ 
4 



1 3Ns 



Hi 



(3.16) 



where u = 1 4- [tanh i In(c^)] and v = tanh|;ln(c^) 
[tanhiln(c2)] . The probability Pp is defined as pp 



Zgy /Z%v where the index Vg indicates that the boundary 



degrees of freedom shown in Fig. 2(a) are fixed 



In the expansion of Z%v and Z^y^ there are two ways 
of taking non-trivial terms. Firstly, similarly to a high- 
temperature expansion of an Ising model's partition func- 
tion, we have terms which are products of squared Ising 

p 

variables r. Moreover, in the expansion of Z^y, terms 
which are products of the Ising variables on the bound- 
ary F (cf. Fig. 2(a) ), are also non-zero, because the Ising 
variables on the boundary arc fixed to a specific value 
(according to the form of the transformation 5, which in 
this case corresponds to the possible different values of 
the Ising variables {t'"}). These terms are classified ac- 
cording to the relative orientation of the appearing r'"" 
variables and contribute differently when the r variables 
are located near one of the four corners of the boundary 
(cf. Fig.[2(i)|). 

We first expand the denominator of the partition func- 
tion, keeping the most relevant term. The second most 
relevant term is of order 0(u^'^'^v'^) ~ 0([lnc^]''): 



Z&v 



C 2 



cosh(- In c^) 



3Ns 



.Ns 



Oil- 



(3.17) 



In expanding the numerator, we have to keep in mind 
that the r variables living on the dual sites next to the 



boundary, shown in Fig. 2(a) are deliberately fixed. So, 
expanding, we have: 



y'^g _ r,Ns-Nr„- 



cosh(- In c^) 
4 



1 3Ns 



.Ns 



1+E 



-(a -ffi) 





(3.18) 



where 



i0 



TT and the variables r are defined on the 
sites of the dual lattice that are pointed out around the 
site of the direct lattice located at r. The path of the con- 
necting line between the sites of the dual lattice denotes 
the type of the term in the expansion. The additional line 
denotes the part of the boundary (cf. Fig. |2(a) I. The 
sum is over all topologically equivalent contributions and 
over all possible locations r along the boundary line. 

What is clear from Eq. (|3.18p is that there are distinct 
contributions in the perturbation series, coming from the 
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corners of the boundary line. Having such a series expan- 
sion, we can derive a perturbation series for the entangle- 
ment entropy, using Eq. ()2.10p . We are going to restrict 
ourselves to terms up to order 0{{v/u)'^) ^ 0((lnc^)*), 
so as to show the existence of a principle in this type 
of expansion, when a correlation length appears. The 
form of the scries expansion depends on the number of 
degrees of freedom on the boundary and the number 
of corners Nc- Here, according to Fig. 2(a) 
consider the physical limit ^ Ny- 
this expansion we find the result 



Nc ^ 4. We 
At low orders in 



Sa 



ln2- 



ln2 















1 




2 





r 

5 / IJ \ 4 
12 \u 



Nc + 



U 



(3.19) 



From Eq. p.l9p it is clear that the correlation length ef- 
fects amount to a renormalization of the term which is 
proportional to the length of the boundary and to a con- 
stant term which scales with the number of corners. On 
the other hand, the topological constant term remains 
invariant. It is apparent that the form of this result is 
a consequence of the structure of the expansion. Hence, 
the form of Eq. (|3.19p will remain unchanged order by 
order in this expansion within its radius of convergence. 
Hence, the value of the topological entropy is a universal 
property of the topological phase. In particular this form 
holds provided the size of the region is large compared 
to the correlation length. This requirement also applies 
to the distance between singular points (corners) on the 
boundary between the regions A and B. 



IV. THE QUANTUM DIMER MODEL ON THE 
TRIANGULAR LATTICE AT THE RK POINT 

In this Section, we show that the global topological 
degeneracy of the topological liquid phase of a quantum 
dimcr model naturally leads to the existence of a topolog- 
ical term in the entanglement entropy. This term is di- 
rectly related to the existence of winding sectors around 
the region of which the entanglement entropy is com- 
puted. 

We are going to focus on the properties of the ground 
state wave function of the triangular lattice quantum 
dimcr model at the RK point^i, even though our argu- 
ments hold more generally. These arguments hold for all 
topological ground state wave functions which have well 
defined loop representations and keep the assumptions 
made in Section |ll] intact. For example, in the general 
class of wave functions studied in Ref . [20, the Temperley- 
Lieb algebra satisfied by the loops ensures that a loop 
cannot be broken into string segments, but can only be 
annihilated. This property, as we will show for the case 
of the dimcr ground state wave function, enforces the 



existence of winding sectors around the region, whose 
entanglement entropy is calculated. 

The dimer model on a triangular lattice at the RK 
point has the following ground-state wave function^: 



(41) 



where C labels configurations of hard-core dimers on the 
plane triangular lattice and the sum is over all configura- 
tions that may be connected through action of the "flip 
term" (cf. Section |ll] and Appendix [A| . This ground- 
state wave function is topologically ordered with a degen- 
eracy that scales as 4^ , Q being the genus of the surface 
on which the model is defined, and can be described by 
a Z2 gauge theory. For example, on a torus, there are 
four equivalent topological sectors, each containing one 
of the four degenerate ground states. In the following, 
we will consider (14.11) defined in one of these sectors. To 





(a) 



(b) 



FIG. 3: (color online) (a) A part of the transition graph be- 
tween a configuration C (light (red)) and a reference columnar 
configuration Co(dark(blue)). It is clear that the transition 
graph is composed typically of loops and isolated links. If 
the configuration C is changed locally, then either the num- 
ber of loops changes, adding or removing isolated links from 
the loops or the relative positions of the loops and the links 
change (without any shape change). (b) Definition of region A 
with respect to region B. The light (green) dimers define the 
configuration in A and the dark(red) ones the configuration 
in B. More generally, any link that lies on the boundary line 
or inside of it belongs to A. 

elucidate more the structure of the ground state wave 
function (j4.1[) . we consider the allowed transition graphs 
of the configurations \C) contained in (|4.ip . with a refer- 
ence configuration |0), say one of the columnar states (cf. 
Fig. 3(a) ). The sum over all possible transition graphs is 
equal to the sum over all possible configurations. Now, 
given one transition graph, we consider all its possible 
transformations, which correspond to all possible gauge 
transformations in the corresponding gauge theory. A 
transition graph typically has two components, loops and 
links, but open lines arc not allowed. The possible basic 
transformations classify to: 

1. Flipping a loop to a set of links, 

2. Shifting a loop with respect to the links, 

3. Combining two loops into a new one or breaking a 
loop into multiple loops. 
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4. Adding a link to a loop or the reverse. 

All of these irreducible transformations correspond to 
flips of different types of flippable dimer loops on the 
actual dimer configurations. Given this construction, it 
is easy to see that there are topological sectors in the 
configuration space of the system. If a loop exists that 
winds around the torus, for example, then no local trans- 
formation out of the four possible ones can lead to the 
removal of this loop, given that this loop is not flipped. 
If the loop is flipped, then the sector has changed. 

Now, we can apply these ideas to the concept of topo- 
logical entropy. Imagine that we divide a closed space 
to regions A and B with a boundary loop F, where A is 
located "inside" B, in the sense that a loop in B which 
winds around A is topologically non-trivial, thus mak- 
ing region B not simply connected. In other words, by 
means of local operations the loop can only be deformed 
to become congruent with the boundary T (cf . Fig. |3(b)[ ) . 
Conversely, a loop in B that does not wind around region 
A is contractible. For clarity, the boundary loop F has 
been chosen in such a way that all dimers that lie on a 
link of the boundary belong to region A and arc defined 
to be degrees of freedom of region A. 

Now, we consider the probability pp, by fixing all the 
degrees of freedom on the boundary F. This can be done 
in this representation by specifying the dimers of region 
B which intersect the boundary loop. Due to the hard- 
core constraint, all possible transformations included in 
the calculation of cannot change the fact that the 
same dimers intersect the boundary at the specified po- 
sitions. We use the construction in terms of transition 
graphs but we choose the reference configuration in such 
a way that all the fixed dimers on the boundary match 
exactly. In this way, we prohibit the possibility that open 
lines are allowed in region B. Then, it is clear that the 
existence of region A introduces a contribution of wind- 
ing sectors of the configuration space (in loop language) 
in the calculation of the probabilities pp of the bound- 
ary configurations. We select a single transition-graph 
loop winding around region A and following the rules de- 
scribed, we can conclude that it cannot be destroyed by 
using the transformation rules 2 or 4, but it could by us- 
ing rules 1 or 3. However, rule 3 can lead to the destruc- 
tion of a non-contractible loop only if it combines with 
another non-contractible loop around region A. Thus, 
if the number of non-contractible loops is even, then all 
these loops can be transformed to contractible ones by 
applying rule 3. On the other hand, if the number is 
odd, then the existence of a single non-contractible loop 
is unavoidable. Then, rule 1 or equivalently, a large trans- 
formation around region A is required to eliminate the 
loop. This fact leads us to the conclusion that 



(4.2) 



where pp and pp represent the probability pp with the 
constraint that there is an even or odd number, respec- 
tively, of non-contractible loops which wind around re- 



gion A. The emergence of this odd/even decomposition 
reveals the topological feature of the underlying gauge 
theory. 

Now, we would like to ask whether there is any special 
relation of pjt with pjT . If the correlation length of the 
system was zero, then the model would become one of the 
loop models considered in Ref. ^2^ with d-isotopy param- 
eter d = 1 (cf. Appendix E]) . In this case, it turns out 
that Pp = pP exactly2£. In the same spirit, for topolog- 
ically ordered ground-states of quantum dimer models, 
systems with finite correlation length, there is no reason 
why configurations in the odd sector should be favored or 
disfavored over those in the even sector, except possibly 
near the boundary F. Therefore, we conjecture that: 



Pp (l + 0(e 



(4.3) 



The reason why Eq. 
the arguments in Ref. 



l|4jP is expected to be true follows 
H, for the degeneracy splitting of 
non-trivial topological sectors in a finite system. One can 
define, in general, large gauge transformation operators 
T, around region A, which commute with the Hamilto- 
nian and any other physical observable operator. These 
operators map a state in one sector (-1- or —) to a state in 
the other. In the dimer language, these operators fiip a 
flippable dimer loop around region A. If we consider the 
operator S which has the property (G| S \G) = pf then 
S" is a sum of products of dimer density operators along 
the boundary contour F. Being in a topological phase, 
S commutes with T, with corrections which depend on 
the tunneling probabilities between the H- and — sectors. 
These probabilities, as shown in Ref. [l|, generally scale as 
e"^r/C^ leading uhimately to Eq.lO 
The quantity pjT can be estimated as follows: Given 
that the length of the boundary is ir, the configuration 
that has the highest probability is the one which con- 
stitutes a flippable dimer loop along the boundary line. 
For a topological ground-state, the probability of such 
a flippable loop is just half the expectation value of the 
Wilson loop on the boundary F and therefore is equal to^^ 
e-'=o(C)ir+ci({)^ where co('f) and ci(^) are correlation- 
length related effects which depend only on the shape 
of the chosen boundary contour F. Therefore: 



,-Co'(?)ir+Ci«(J) 



(4.4) 



The argument holds for any pp so thus the entanglement 
entropy of region A is: 



Sa = -^2pf^ln(2pf^) 

9 

= cO(C)ir-ci(0-ln2 



(4.5) 



where c0^i(0 = E,Pr,c;^;'(C)- 

This type of argument can be repeated for any similar 
topologically ordered quantum dimer model ground state 
wave function like the ones on the Kagomoiii and Fisher— 



11 



lattices. The existence of large transformations winding 
around region A introduces non-trivial winding sectors in 
the configuration space with a non-trivial contribution to 
the calculation of the entanglement entropy. 



V. DEFECTS, TOPOLOGICAL DEGENERACY, 
AND TOPOLOGICAL ENTROPY 

Defects, or equivalently violations of the local con- 
straints of the model, can be added in both dimer and 
vertex models in a similar way, given that there are 
known mappings between these classes. Such defects can 
be regarded as representing fluctuations of a matter field 
which, in the cases at hand, carries a Z2 charge. We are 
going to focus on defects added in the ground-state wave 
function of a deconfined topological phase at the Kitaev 
point in two different ways. Firstly, we examine the effect 
of such defects on the eight-vertex quantum model at the 
Kitaev point. In this model, in the representation where 
the link operator "electric field" is diagonal, (electric) 
defects are vertices which have an imbalance of ingoing 
and outgoing arrows, namely vertices with 3 arrows in, 1 
out or 3 out, 1 in. (See the notation and terminology of 
Section Hill and Appendix 1X1 ) On the other hand, in the 
representation where the link operator ( "gauge field" ) f'^ 
is diagonal, (magnetic) defects are defined by plaquettes 
which cannot resonate. Next, we examine the effects of 
virtual magnetic and electric charges on the Kitaev state 
by looking at these effects in the problem of a Z2 gauge 
theory with matter. These two cases behave in clearly 
different ways which illuminate the problem of the sta- 
bility of the topological phase. 

It is clear that there are two physically distinct ways 
of adding such defects in the ground state wave func- 
tion. Firstly, one can allow defects to be mobile, like a 
dilute gas. In this case, even though the ground-state 
wave function is connected perturbatively to the Kitaev 
wave function, local but non-perturbative terms have to 
be included in the Hamiltonian of the quantum model^. 
Such a wave function corresponds to a state in which 
the Z2 matter field has "condensed" , in the sense that 
free (Z2) gauge charges are proliferating in the ground 
state. It is a well known result^ that the state in which 
the matter field condenses is smoothly connected to a 
confining state. In such a state, the loop configurations 
are essentially broken up into strings of finite length (the 
confinement scale). In such a state, the loops cannot ex- 
plore the topology of the manifold on which the system is 
defined, or of the region of space being observed. Hence, 
such a state does not represent a topological phase. 

On the other hand, the quantum Hamiltonian can be 
perturbed by a term that favors the existence of such 
defects and charge neutrality should enforce the emer- 
gent defects (either magnetic or electric, see Appendix E| 
to be created in pairs. In perturbation theory, the cor- 
rections to the ground state wave function contain even 
numbers of defects which are separated mutually by a 



distance typically of the order of the term in the expan- 
sion series. Only at very high orders in the perturbation 
series are defects more or less mobile, and an interpola- 
tion between the two limits of free and confined defects 
could be possible, if a phase transition is not present in 
between. Given that such a transition is typical in Ising 
gauge theories with matter fields in 2 -I- 1 dimensions^, 
one should expect that a similar transition happens in 
this case as well. 



A. Defect liquids 

We will now consider an extension of the Baxter wave 
function discussed in Section IIIII to include defects that 
violate the eight-vertex constraint. Namely, we will allow 
configurations with three arrows in and one out and vice 
versa. We will follow the notation and terminology of 
Section IIIII and Appendix [A] 

More specifically, we will write the Kitaev wave func- 
tion as a limit of a state which allows for the constraint to 
be violated at every vertex {r} of the lattice. In particu- 
lar, if the Hilbert space is enlarged to include all possible 
arrow configurations, then the Kitaev wave function can 
be written as follows: 

\G) = ^e^^^-^"'"^-"""'-!) |{t}) (5.1) 

{r} 

where i,j,k,l label the four links surrounding the ver- 
tex r, and we have denoted by t; the eigenvalue of the 
electric field operator (= of Section III) on link i 
(and similarly with j, fc, I). By direct examination of the 
wave function of Eq. (|5.ip it is straightforward to see 
that in the limit K-^ 00, the exponential term enforces 
the eight-vertex constraint (even number of arrows in- 
even out) on each vertex of the lattice by setting to zero 
the amplitudes of all Ising (arrow) configurations which 
violate the constraint. It is instructive to generalize a 
little more this wave function by adding an analogue of 
the classical polarization field in the classical eight-vertex 
model: 

\G) = ^e^.-f^"(n^.^^^<-i)+^*'E.n |{^}^) (5 2) 

{r} 

When Kg 1, the only configurational state that has an 
appreciable amplitude is the one where = -|-1 on every 
link i of the lattice, a "ferroelectric" configuration. On 
the other hand, when Ki 0, the usual Kitaev point 
is approached where the state is in a liquid phase. It is 
clear that a phase transition should separate these two 
limits and it does. 

It is immediate to see that the norm of the state \G) 
of Eq. (|5.2p is simply the partition function of the 1 + 1- 
dimensional Ising gauge theory coupled to a matter field. 
A well known duality transformation for Ising gauge the- 
ories in 1 + 1 dimensions^i^, which maps the 1 + 1 dimen- 
sional (in an Euclidean space-time lattice) Ising gauge 
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theory with a matter field to a classical two-dimensional 
Ising model with an external magnetic field, can be ap- 
plied here. Under this duality transformation, the norm 
of the ground state (|5.2p is written in terms of Ising vari- 
ables Si, defined on the vertices of the lattice (instead of 
the links). Following closely Refs. [H and[3l|, one easily 
finds that the norm of the ground state is: 



(G|G) 



1 



\2N, 



(cosh(2if„))^=(cosh(i^f)) 

^e^<'^>'^*"'''^+'''^'-'" (5.3) 

{s} 



where i and j now denote sites of the dual lattice, and 
/?* = - i In tanh 2Ke, h* = -^\n tanh 2A',, (5.4) 



It is important to notice that the factor 1/4 in Eg. ([57 
comes from the fact that on a torus, where the system is 
defined, periodic boundary conditions should be imposed 
on the direct lattice. This means that on the dual lattice, 
products of pairs of dual variables should equal on the 
boundary, leading to a double degeneracy for each of the 
directions of the torus. 

The amplitudes of the states that compose the ground 
state wave function follow the distribution of an Ising 
model in a magnetic field. In the defect-free limit Ky — > 
oo, it follows from Eq. that h* = 0. If h* ^ 0, then 
defects are present and independent from each other. So, 
the presence of mobile defects in the ground state of the 
quantum eight vertex model corresponds to a magnetic 
field which breaks the Ising symmetry of the weights of 
the wave function. On the other hand, if we also set Kg ~ 
(which implies /3* oo) to reach the Kitaev point, 
the weights of this wave function now have the form of 
the Gibbs weights of an Ising model at zero temperature 
and magnetic field. In this limit, there are just two dual 
configurations which contribute a non-zero amplitude to 
the ground-state structure, the two fully ferromagnetic, 
up and down, Ising configurations. It is clear that this 
degeneracy in this dual formulation of the ground state 
wave function corresponds to the topological order that 
characterizes this wave function. 

Let's consider the calculation of each as it is de- 
fined in Section Hill and Fig 2(a) i.e. the probability of 



a configuration on a loop F, the boundary of the region 
A under observation for the calculation of the entangle- 
ment entropy. Using the methods of Ref. [33l . an Ising 
variable in a correlation function maps under duality 
to the operator e~'^^ where (i*,j*) represent the 

sites whose midpoint is the site i. The expectation value 
of such an operator in the dual formulation, in the limit 
/?* — > oo, /i* = 0, is clearly zero, because antiferromag- 
netic configurations are not present. In this way, 

where the distribution of ±'s in Eq. (|5.5p . depends 
on the choice of the loop configuration gi-. In the 




(b) 



FIG. 4: (a) The positions of the arrows denote the direct lat- 
tice and the dots (located on the vertices) the corresponding 
dual lattice. The dashed line marks the boundary of region 
A and the thick bonds are the bonds on which the dual cou- 
phng (3* has to reverse its sign when a correlation function 
(fir), a dual Wilson loop, is to be computed, (b) The string 
representation of a typical eight-vertex configuration. 



relevant limit /3* oo, all but two terms of the 
numerator's expansion are zero. The only non-zero 
terms are the first one 2"^"" (G |G) and the last one 
2-Nr ^(;;| j-|.^^g-2/3*s..s,. ]^qy rj,^^ Ya.ttei term is non- 
zero and equal; in absolute value, to the former one be- 
cause of the following reason: when all bonds on the 
(dual) loop F are antiferromagnetic, they can all be sat- 
isfied without any energy penalty by having all spins in 
region A being -1-1 and all spins outside region A being 
— 1, or the reverse. If the term is negative, then pp = 0, 
as it is expected because whenever Y\tx = —1 around 
the loop, the eight-vertex constraint is violated. So, in 
the hmit /3* = oo, ft,* = we have: 



g _ r,l-Wr 



(5.6) 



for all the non-zero probabilities, and the known result 
in Eq. p.l3|) is recovered, as expected, leading to the en- 
tanglement entropy in Eq. p.lSp . 

If /3* 3> 1 but finite(/i* = 0), then it is evident that the 
degeneracy between the two dual ferromagnetic states is 
present and exact until /3* = /3c = ^ ln(V2 + 1), where an 
Ising transition takes place to a paramagnetic state. For 
[3* > Pc, the topological entropy, using the same argu- 
ments as above, remains constant and equal to In 2 until 
/?* = (3c where it drops discontinuously to zero, since 
the exact two-fold degeneracy, which led to the emer- 
gence of the topological entropy, is not present (see also 
Ref. HI)- For (3* > Pc, the constant term of the en- 
tanglement entropy, computed from a low-temperature 
expansion of the dual system for /?* ^ 1, contains expo- 
nentially small corrections(~ e"*^^ ), coming from anti- 
ferromagnetic fluctuations near the corners of the bound- 
ary between the ferromagnetic regions A and i?(cf. Sec- 
tion |lll|. At the other limit, also similarly to our results 
in Section IIIH for f3* 0, only constant terms coming 
from the corners between the disordered paramagnetic 
regions A, B are present, not related to the topology of 
the system. 

Given this simple way of calculating correlation func- 
tions, we can ask what is the effect of Ky ^ 1 but finite. 
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or equivalcntly h* ^ Q, introducing mobile defects in the 
wave function. Given that the exact degeneracy of the 
two dual ferromagnetic states is lifted when /i* 0, it is 
clear that the topological term should not be present in 
this case. This is actually the case in the limit /?*—!■ oo 



q 

Pr 



^ cosh[{NA - NB)h* 



sh{Nsh*) 



(5.7) 



where Na and Nb are the numbers of dual sites in regions 
A and B respectively. The entanglement entropy is then, 



Sa = {\n2)Nr + 



/cosh[(A^A - Nb W 
\ cosh{Nsh*) 



(5.8) 



In the thermodynamic limit, clearly Sa = (ln2)7Vr, and 
the topological entropy has clearly been wiped out by the 
existence of mobile defects, i.e. -f — 0. 

On the other hand, when /?* ^ 1 but finite, clearly 
effects similar to the ones described in Section |TTT] natu- 
rally appear. The reason is that on the direct lattice, the 
term proportional to Ki in the amplitude weight of the 
wave function corresponds to a "ferroelectric" field in the 
classical eight vertex field, which, when weak, has simi- 
lar effects to the ones described in Section |lll](c.f. classic 
results from Baxter—.) In terms of the dual formulation 
we presented here, it is clear that the exact Ising degen- 
eracy, which leads to the topological entropy, persists for 
h* — until (3* ~ /3*, where an Ising transition to a 
disordered state takes place. 



B. Perturbation theory and lifting of the 
topological degeneracy 

The case <g; 1 , —f oo can be viewed also in 
another way. One can perturb the Kitaev Hamiltonian 
Eq. (|3.3p with the following perturbation: 



(5.9) 



where < X <^ t. The effects of such perturbations 
have been studied recently, using numerical methodsi^i^ 
who found, as expected, that the topological (deconfined) 
phase is stable. As it is discussed in detail in Appendix 
lAl this operator corresponds to an electric field in the 
language of the Ising gauge theory. If we consider the 
string representation of the eight-vertex configurations 
(cf. Fig. [4(b)]), which can be defined by putting strings 
on links with an arrow pointing up or right, then clearly 
the term (|5.9p does not favor the existence of strings on 
links of the lattice. In this representation, similarly to the 
discussion in Section IIV| topological sectors on a closed 
topological manifold, such as a torus, are distinguished 
by whether there is an even or odd number of strings 
winding each of the directions of the torus. Due to this 
fact, the perturbation Eq. (|5.9p hfts the topological de- 
generacy, making the (-I-, -I-) sector the most favorable, 
where the state with no strings belongs. 



From another point of view, the operator — 1 cre- 
ates two magnetic defects with an energy penalty 2t in 
the ground-state subspace, if one views the ground-state 
wave function in the representation where is diagonal. 
We will consider the states created by the application of 
products of these operators on the ground-state as the 
basis for applying perturbation theory. By applying for- 
mal Brillouin-Wigner perturbation theory, we take: 



|G) = \Gk) 



E 



G 



{r} 



V 



G 



G 



E- 



\G 



K 



{■■} 

ci IGi) 



G 



where |Gi) = Y.{v} 



G 



(1) 

{r} 



and ci = 



2t+A ■ 



(5.10) 

(5.11) 

The state 



G|^j ) denotes the superposition of all possible configu- 
rations with 1 magnetic defect at the location r. Also, it 
is clear from the structure of the perturbation series that 
at a higher order n, the correction |G„) emerges with a 
coefficient c„ of the order 0((A/i)"). 

The entanglement entropy of the above wave function 
can be computed through our usual procedure, but only 
terms up to order 0[c\) are going to be consistently cor- 
rect. Higher orders in the entropy require the knowledge 
of higher orders in the perturbation series. Given that 
{Gk\ |Gi) ~ 0, we may consider the perturbed ground 
state as a pure state with unit amplitude (unnormalized) 
for 0-defect configurational states and ci for one-defect 
ones. Then, one can define a partition function (as in 
Section HIT)) and associated probabilities Pp, and proceed 
to the calculation of the entanglement entropy, as in Sec- 
tion [nil Following similar steps as before we have, 

Sa = {\n2 + cl + 0{{\/tt))NT-\n2 (5.12) 

It is clear that the topological term is robust under such 
low order perturbations and also, higher orders lead to 
corner effects, in a similar way as it is discussed in Sec- 
This behavior has been verified numerically in 

From the form of the perturbation Eq. (|5.9p . it is clear 
that the exact topological degeneracy of the winding sec- 
tors, discussed in Sect ion [TVl and related to the string 
representation (cf. Fig |4(b)[ ), is lifted. This happens be- 
cause, if one considers the perturbation applied only on 
a loop around region A, then in the -|- sector, the state 
with zero strings winding around region A is favored, but 
in the — sector, the state with 1 string is favored with 
different energies. 

Given that the degeneracy is lifted, it is clear from 
our discussion throughout this paper that the entangle- 
ment entropy should acquire a non-trivial correction to 
the constant topological term, not related to the shape 
and type of the boundary, but directly related to the lift- 
ing of the winding sectors' degeneracy. In our example, 
this correction is identifiable. In order to show the prin- 
ciple, let's keep just one relevant additional term in the 
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wave function, which amounts to a set of terms which 
emerge in the perturbation series of Eq. (j5.10p to high 
orders: 

g) = |Go) + c„^, n(^'T' - 1) |Go> (5.13) 
r' 

where Cn^, is a constant proportional to (A/t)"r" where 
nr' the length of a generic chosen loop which winds re- 
gion A. The correction to the original wave function has 
the property that is non-zero whenever a string occupies 
all the links which lie on the loop F'. On the other hand, 
in the representation where is diagonal, the correction 
should be interpreted as the creation of a line of magnetic 
defects. Given that the winding sectors of the configu- 
rations surrounding region A are characterized by even 
or odd number of strings, it is clear that this correction 
to the wave function has topological features. We can 
apply similar techniques as before to calculate the entan- 
glement entropy, given there is no complexity with the 
definition of the boundary degrees of freedom (the loop 
r' is chosen to be generically away from the boundary T 
of region A) . We can calculate the entanglement entropy, 
leading to the following results: 

2^-i + (2c„^, +c2^,)2^-^r' ^^-^^^ 
= 2i-^-(l - (2c„,, -I- <,)2-^r- + . . . ) (5.15) 

and the entanglement entropy is: 

Sa = (In 2)Nr - In 2 + 2-^r' (2c„^, + c^^, ) + ■■■ 

(5.16) 

From the above expression it is clear that a correction 
to the topological term 7 = ln2 appears, of the order 
~ (A/(2i))^r' . This is an exponentially small correction 
to the topological entropy which vanishes (exponentially 
fast) in the limit of a large region A, provided the per- 
turbation theory used here is convergent. This is indeed 
correct for A < Ac, which is to say within the topological 
(deconfined) phase. On the other hand, these topological 
excitations, when they proliferate (A ~ Ac), ultimately 
lead to the destruction of the topological entropy and, 
for a homogeneous system, to a global topological phase 
transition. As it is well known, the physics of this quan- 
tum critical point at Ac is equivalent (upon duality) to 
that of the quantum critical point of the 2-1-1 dimen- 
sional Ising model in a transverse field, a.k.a. the clas- 
sical three-dimensional Ising model. The quantity of in- 
terest here, the topological entropy, is thus related to the 
statistics of the proliferating domain walls of the 3D ising 
model at its critical point. The solution of this problem 
is still open. 

Finally, we note that in this Subsection we considered 
only the effects of fluctuations on magnetic charges. The 
methods used here can be applied to the case of electric 
charges as well. Moreover, the well known self-duality 



of the Ising gauge theory with Ising matter in 2 -|- 1 
dimensions^i^, discussed in Appendix [3 implies that 
this result also applies for electric charges, again within 
the domain of convergence of perturbation theory about 
the Kitaev (deconfined) point. This perturbation theory 
is well known to have a finite radius of convergence. We 
thus conclude that the topological entropy 7 = In 2 is a 
property of the entire deconfined phase and not just of 
the Kitaev limit. 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

The calculations presented in this paper have implica- 
tions for numerical attempts to determine the topological 
order in a ground-state wave function in the case of a fi- 
nite correlation length. The first issue to note is that 
there will be non-universal corrections to the sublcading 
piece, as noted in Eq. (|2.12p . and in the general case, the 
function 6(^) will not be known. This makes the prob- 
lem of extracting 7 more subtle than just calculating Sa 
versus L and then plotting the intercept. 

However, the observation in Section Hill that such non- 
universal corrections come from the corners of region A 
suggests a way around this issue. For example, on a 
square lattice and assuming the definitions of the bound- 
ary and its degrees of freedom of Section IIIIl the en- 
tanglement entropy of a region with four corners will 
have the scaling form = a4(^)i -I- 64(^) -I- ... where 
64 (^) = 4c(^) — 7 and c(^) is the contribution of a sin- 
gle corner which will depend only on whether the corner 
lies on the direct or dual lattice. The entropy S*^^ of 
a region with six corners will have a similar form where 
66(<^) = 6c(^) — 7; here we have used the fact that a 90° 
and 270° corner give the same contribution. Therefore, 
if we calculate S^"* and S^"* versus L and determine the 
respective intercepts, then 7 may be extracted by: 

7 = 266(0-364(0 (6.1) 

This type of subtraction to eliminate the unknown effects 
of corners was one of the factors motivating the construc- 
tion of Ref.0. The present example is, perhaps, a simpler 
example of the same strategy. It is clear, though, that 
the application of such methods requires a precise under- 
standing of the relevant degrees of freedom and the pre- 
cise definition of the boundary shape and type between 
regions A and B. 

Additional simplifications arise for the special class of 
wave functions emphasized in this paper. In Section [iTl 
it was shown that computing the entanglement entropies 
of wave functions whose normalizations resemble parti- 
tion functions of classical statistical mechanical systems, 
amounts to computing probabilities of boundary config- 
urations of the same classical systems. Such probabilistic 
quantities are convenient to compute using classical ap- 
proximate techniques such as Monte Carlo simulation. 
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The advantage of such techniques is that it becomes pos- 
sible to investigate system sizes large enough that in- 
trinsic finite size effects, such as corners interacting with 
each other or perturbative mixing of topological sectors 
as discussed in Section IV Bl are no longer an issue. 

Moreover, by working with the reduced quantities men- 
tioned in Eqs. (|2.13|1 and (|2.14p . the problem further re- 
duces to one of computing expectation values of non-local 
loop operators of varying loop shapes. The topological 
entropy may be extracted using the strategy of Eq. (j6.ip . 
for example, though care should be taken to choose the 
loop operators such that the different corners of the loop 
have similar immediate environments so that the sub- 
tractions go through. 

This viewpoint provides a route to answering certain 
unresolved questions about Fig. 1(b) away from the Ki- 



taev point such as the mechanism by which the topologi- 
cal order is lost along the critical lines to give way to the 
ordered "antiferroelectric" phases. Such a study, which 
further connects the present work to ideas discussed in 
Ref. [H, is a natural topic for further investigation. 

In this paper, we tried to elucidate the topological 
features of the entanglement entropy of topologically or- 
dered ground-state wave functions with finite correlation 
length. Firstly, by expanding generally the entanglement 
entropy in terms of appropriately defined boundary prob- 
abilities, we suggested that reduced quantities, directly 
related to the usual von Neumann entanglement entropy, 
contain the same topological information as the entangle- 
ment entropy. Then, we showed explicitly in the context 
of the quantum eight- vertex model, that the effects from 
the existence of a finite correlation length amount to con- 
tributions to the linear part and non-topological constant 
pieces coming from non-smooth parts of the boundary. 
Then, with reference to the topologically ordered ground- 
state wave function of the quantum dimer model at the 
RK point, we showed that the concept of topological en- 
tropy is connected to the existence of ficticious sectors 
in the calculation of the entanglement entropy, allow- 
ing for the identification of the topological order of the 
state. Finally, we considered the effects of topological de- 
fects on the entanglement entropy of the quantum eight 
vertex model near the Kitaev point. We firstly consid- 
ered the non-perturbative regime of a defect liquid, where 
defects are mobile and cause the topological entropy to 
vanish. Secondly, we showed that if defects are slightly 
favored in the quantum Hamiltonian. the topological en- 
tropy is a topological invariant and effects related to the 
shape and type of the boundary are only allowed. Finally, 
along these lines, we identified the terms in the pertur- 
bation expansion which render topological corrections to 
the topological entropy and should ultimately lead to a 
global topological phase transition, if they proliferated. 

Note: As this work was being completed, we became 
aware of the work by Castelnovo and Chamon^i^ who in- 
dependently, and among other questions, considered the 
wavefunction in Eq. (|5.2p in the limit — > oo and stud- 
ied its topological entropy as a function of in the con- 



text of an investigation of the stability of a topological 
phase. We thank these authors for communicating their 
results with us before publication. 
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APPENDIX A: A SUMMARY OF RELEVANT Za 
GAUGE THEORY RESULTS 



In this Appendix we will review a number of stan- 
dard and well known results from the gauge theory lit- 
erature in the context of its applications to topological 
phases. Much of what we discuss here was developed 
and reviewed extensively in Refs.l36l.l37l.l38l andl30l. The 
topological nature of deconfined phases of discrete gauge 
theories was emphasized in Refs. [S^ and and its role 
in topological mechanisms for quantum computing was 
first formulated by Kitae\*^ (for a recent review see l4ll.) 
Here we will focus on the deconfined phase of Z2 gauge 
theory with matter— i^i^i^l which is the simplest exam- 
ple of a topological phase. Some aspects of this problem 
were revisited recently in Ref. [4^ . 

Let us consider a Z2 gauge theory with matter fields 
on a 2 + 1-dimcnsional square lattice in the Hamiltonian 
formulation. As usual the Z2 Ising gauge fields live on 
the links of the square lattice. We denote by the Pauli 
matrix T|(r), the gauge field on the link (r, r + ej) of 
the square lattice (with j = 1, 2 denoting the two spacial 
directions). The (Ising) matter field resides on the sites 
r of the square lattice and is denoted by the Pauli matrix 
^^(r). The quantum Hamiltonian is 



H = -g ri(r)-^c7.(r) 



rj = l,2 



-^ri(r)rf(r + ei)ri(r + e2)r|(r) 

r 

-A ^ (T^(r)Ti(r)cr^(r + ej) 

rj = l,2 



(Al) 



where g and A are two coupling constants. The first term 
in Eq. (jAl[) is a gauge field kinetic energy (electric field- 
like) term, the second is a matter kinetic energy (a trans- 
verse field), the third term acts on plaquettes and is a 
gauge potential energy (magnetic flux) term, and finally 
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the last term minimally couples the matter and gauge 
fields. 

This quantum Hamiltonian has a local Z2 gauge sym- 
metry. The generators G(r) of local Z2 gauge transfor- 
mations are, 

G(r) EEa.(r)ri(r)ri(r-ei)r,2(r)r,^(r-e2) (A2) 

These are Ising operators 

G(r)2 = 1 (A3) 

which commute with each other 

[G(r),G(r')]=0 (A4) 

and with the Hamiltonian 

[G(r),iJ]=0 (A5) 

The physical Hilbert space are the gauge invariant states 
iPhys), 

G(r) IPhys) = IPhys) , V r (A6) 

i.e. the physical states invariant under arbitrary local 
time-independent gauge transformations. In the basis in 
which the gauge field r| (r) and the matter field cr^ (r) are 
diagonal, the action of the local gauge transformations 
G(r) is to flip the sign of the matter field at site r and 
of the gauge fields on the four surrounding links. The 
simpler physical states of this theory can be qualitatively 
described as magnetic charges, i.e. plaquettes where the 
magnetic flux term is —1, and electric charges, i.e. sites 
r where crx(r) = 

A local gauge symmetry cannot be spontaneously 
broken.— However it is possible to use the gauge invari- 
ance of the theory to fix the gauge. In this case it is 
possible to use the action of the gauge generators G(r) 
to fix the gauge, 

a.(r) = 1 (A7) 

globally and completely, on all sites r of the square lattice 
for any boundary conditions (a disk, a torus, etc.) This 
is the unitary-London gauge. In a general gauge theory, 
the gauge in which the phase of the matter field is set to 
zero is known as the unitary gauge, and in the theory of 
superconductivity it is known as the London gauge. In 
this Ising theory, this gauge fixes the matter degrees of 
freedom completely. 

In this gauge the Hamiltonian takes the form 

rj = l,2 

- ^xi^Hir - ei)T2(r)rJ(r - 62) 

r 

-Eri(r)r|(r + ei)ri(r + e2)r,2(r) 

r 

-A J2 

r,i=l,2 



In this gauge the local symmetry is absent (which is nat- 
ural since we have fixed the gauge completely). 

Ising gauge theories, with and without matter, have 
well known duality mappings (^i^i^i which arc straight- 
forward generalizations of the Kramers- Wannicr duality 
of the 2D classical Ising model. In short, in 1 + 1 dimen- 
sions Ising models are self-dual, while in 3 + 1 dimensions 
Ising gauge theories are self dual instead. In 2 + 1 di- 
mensions Ising models are dual to gauge theories, and, 
in particular, quantum Ising gauge theories coupled to 
Ising matter fields are self dual in 2 -I- 1 dimensions. In 
the context of the unitary-gauge Hamiltonian, Eq. ljASp . 
self duality of the theory amounts to the simple mapping 
g X, which requires the phase diagram to be invariant 
under this symmetry, with the caveat that the electric 
and magnetic properties of the states are exchanged un- 
der this duality symmetry. 

The phase diagram of this theory is well knowni^ If 
the dimensionality d of space is greater than 1 , it has two 
distinct phases: a) a confined phase, smoothly connected 
from a so-called Higgs phase (from which it cannot be 
distinguished) , and b) a free charge (or Coulomb) phase. 
The confinement-Higgs phase has been extensively dis- 
cussed in the literature^ and we will not discuss it here. 
It is a massive phase with a unique ground state and 
a gauge-invariant spectrum of massive magnetic charges 
and no electric charges. It occupies the strong coupling 
sector of the phase diagram, where the gauge coupling g 
and/or the matter coupling A are large. 

The deconfined or free charge (Coulomb) phase occurs 
when both the gauge and matter couplings are weak, 
< .9 < 5c and < A < Ac, where gc and Ac are 
critical couplings. In this phase the spectrum contains 
finite energy states which carry the Z2 charge^. The de- 
confined phase is controlled by the infrared stable fixed 
point ai g = X = which describes a gauge theory in the 
extreme deconfined limit and matter fields which are in- 
finitely heavy. At this fixed point the Hamiltonian takes 
the simpler form 

H = -Eri(r)ri(r-ei)r2(r)r2(r-e2) 

r 

- ri(r)r|(r + ei)ri(r + e,)rUr) (A9) 

r 

known as Kitaev's toric code^ Hamiltonian. It consists 
of the sum of two mutually commuting sets of operators: 
the first term is the generator of time-independent gauge 
transformations of the pure gauge theory and it is the 
product of four operators residing on the links em- 
anating from the site (or vertex) r. The second term, 
in which the operators Tz reside on the links of the pla- 
qucttc labeled by the site r, is the actual Hamiltonian. 
Clearly, in this limit, the eigenstates of the first term 
are states with or without Z2 electric charges (present if 
this operator takes the value —1) while the eigenstates of 
the second (plaquctte) term arc states with or without 
Z2 magnetic charges (present if this term takes the value 
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— 1). In this extreme deconfined limit, the ground state 
has neither electric nor magnetic charges. 

We end this Appendix with a description of two alter- 
native, dual, descriptions of the states in the deconfined 
phase. 

1. The first approach is the standard description of the 
states in a gauge theory in terms of the eigenstates 
of the vector potential which in this Z2 theory are 
the eigenstates of the operators T^ (r), which are 
not gauge-invariant states. At this stage there are 
two options on how to proceed. One option, which 
is indeed the standard approach, consists of fixing 
the gauge. This works provided it is possible to fix 
the gauge completely, something which is not pos- 
sible in many gauges on manifolds with non trivial 
topology such as the torus. The alternative is to 
define gauge invariant states as the linear superpo- 
sition obtained by the action of the gauge group 
Z2 on some reference gauge-fixed state, i.e. by all 
possible actions with the gauge generators G(r). 
In this approach it is essential to account for the 
action of the large gauge transformations, transfor- 
mations in which the action of the gauge generators 
wraps around the non contractible closed contours 
Ffe of the space manifold, e.g. k = 1,2 for the 
torus. This effectively splits the configurations of 
gauge fields into classes each labeled by the expec- 
tation value of the (gauge-invariant) Wilson loop 
operators along each non contractible contour Tj, 
i.e. Wj = rirer '^z('')- ^he other hand, the 
generators of large gauge transformations are the 
operators that create Dirac strings along the same 
non contractible contours Tj, the closed contours 
on the dual lattice. These Dirac string, dual, gen- 
erators are denoted by Wj = Orer These 
mutually dual sets of Wilson loops obey the simple 



algebra. 
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), W^fcj = 0, for j ^ k (AlO) 



2. The second alternative consists in working directly 
instead with gauge-invariant states, as discussed 
extensively by Freedman and coworkers This is 
the approach that we followed in this paper. In 
this approach, the space of states are eigenstates 
of the Wilson loop operators on arbitrary closed 
contours F, contractible or not. In addition, one 
uses the eigenstates of the electric field link opera- 
tors r^, which are gauge invariant by construction. 
A simple picture of the space of states in this ba- 
sis is obtained by mapping it to the set of closed 
loop configurations on the square lattice: a link be- 
longs to a loop if Tj; = — 1 on that link and it does 
not otherwise. The magnetic plaquette term of the 
Hamiltonian acts on neighboring loops and recon- 
nects their strands (thus changing the loop config- 
urations). When acting on an elementary closed 
loop configuration, which occupies just one plaque- 
tte, it annihilates it with weight d ~ 1. This alge- 
bra furnishes a representation of the generator of a 
Temperley-Lieb algebra with d-isotopy parameter 
d ~ 1, as discussed in Ref. H^. In this picture, the 
ground state in the extreme deconfined (Kitaev) 
limit is an eigenstate of the magnetic plaquette 
term and it is thus the state with equal amplitude 
superposition of all loop configurations, discussed 
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